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Abstract 

Spatial intensity correlations between waves transmitted through random media are analyzed 
within the framework of the random matrix theory of transport. Assuming that the statistical 
distribution of transfer matrices is isotropic, we found that the spatial correlation function can be 
expressed as the sum of three terms, with distinctive spatial dependences. This result coincides 
with the one obtained in the diffusive regime from perturbative calculations, but holds all the 
way from quasi-ballistic transport to localization. While correlations are positive in the diffusive 
regime, we predict a transition to negative correlations as the length of the system decreases. 

PACS numbers: 42.25.Dd,05.40.-a,72.15.Rn 



1 



When a wave propagates coherently through a random medium important correlations 
emerge between the different propagating paths, which manifest themselves as correlations 
in the intensity speckle pattern. They have been the subject of great interest over the last 
decade for the case of temporal, angular, and frequency correlations. l[ 0, 0, U, 0, H> E , H, Bl - 
Recently, the direct observation of spatial correlations in the intensity speckle pattern 10Ull| 

waves transmitted through a random medium 

has renewed the interest in this problem [l^ . 

One of the theoretical approaches followed to study this problem involved a microscopic 
diagrammatic calculation 2> Q 0|. This was also the approach in the recent work that 
showed that the spatial correlation function of the normalized intensity can be expressed 
as the sum of three terms, which differ in their spatial dependence Il0|, and in the work 

n 

finding an equivalent structure for the correlations in the polarized radiation [12J. While 
these diagrammatic approaches have the appealing advantage of illuminating explicitly the 
nature of the correlations in terms of the microscopic trajectories of the underlying paths, 
their application is strictly limited to the diffusive regime. 

An alternative approach, macroscopic in nature, has been applied successfully to study 
angular correlations in random media □ □ □ n ; it considers the correlations between 
the transport coefficients in the scatte ring matrix describing the system in the framework of 
Random Matrix Theory (RMT) [tIQjQ]. Most of the work based on RMT has been focused 
on the study of angular or channel-channel correlations. It shows that only the assumption 
of isotropy of the transfer matrix distribution, discussed below, determines the structure of 
the correlations as a function of channel indices. It is the purpose of this work to apply the 
RMT approach to study the spatial intensity correlation functions. Just from the isotropy 
condition we derive the structure of the spatial correlations of the normalized intensity. Our 
result, which coincides with the one obtained from microscopic perturbative calculations 
(in the diffusive regime) is not perturbative and holds all the way from quasi-ballistic 
transport to localization. Only the specific values of the three coefficients (C\, C2 and C3) 
depend on the transport regime. These values can be obtained from the Monte Carlo solution 
of the Dorokhov, Mello, Pereyra, and Kumar (DMPK)Q, Q] scaling equation, and are in 
full agreement with microscopic numerical calculations of bulk disordered wires. While the 
correlations are positive in the diffusive regime, we predict a transition to negative values 
for both angular and spatial correlations as the length of the system decreases. 
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We will consider a wave propagating in the ^-direction in a constrained geometry (see 
the inset of Figure 1). The input and output faces of the disordered region are the planes 
z = and z = L respectively. The transverse coordinates in the system are described 
by ~p. We will only discuss the case of scalar waves in this work, neglecting polarization 
effects. The eigenfunctions of the cavity (in the absence of disorder) naturally separate into 
a longitudinal and a transverse part, 

0n(j) = -=^n(p) exp {±ik n z} (1) 

V fen 



The integer n = 1,2, ... ,N labels the propagating modes, also referred to as scattering 




channels. Mode n has a real wave number k n = \Jk 2 — q 2 , , where k is the wavenumber of 
the incident radiation, and q n is the momentum associated with the normalized transverse 
wave function ?/> n (p). The normalization of the total wave function <p n is chosen to carry 
unit current. The scattering matrix S 

(2) 

relates the asymptotic propagating outgoing waves to the incoming ones: Matrix elements 
Tf, a and tj a denote the reflected amplitude in channel b and the transmitted amplitude in 
channel j when there is a unit flux incident from the left in channel a; and have 
an analogous meaning when the incident flux in channel i comes from the right. Flux 
conservation and reciprocity imply that the matrix S is unitary and symmetric. 

Let us consider two point sources at ~r = f\ and 
~r = fs on the left hand side of the system, and two detectors at ~r = r[ and ~r = 7% on the 
right side (as sketched in Fig. 1). For a point source at ~r = f\ the incoming field from the 
left is proportional to the Green function of the clean waveguide, 

Gt^A,r)= l -Y^ti\r A )<P + a (r) ; (z > z A ), (3) 

a 

The field at a point Ti on the right side outside the system will be given by 



(4) 



where c a oc 4>+*(j a)- The average intensity at that point is given by {I (A, 1)) = (\E(A, 1)| 2 ), 
where < ... > denotes disorder averaging (over the ensemble of samples). For a finite sample, 



the transmission amplitudes, tj a , are assumed to have random phases with (tj a ) = and 
(tjat*t a ,) = (Tj a ) 5 aa i5jji . The square of the field-field spatial correlation function can be 
written as 

cl(A1,B,2) " (WW (5) 

where 

(E(A,1)E*(B,2)) = Y,Ca<'rt(nW(^){t ja t* fa ,) 

aa'jj' 

= J2 c Mj(ri)<t>r(r2)(T ja ). (6) 

(with d a oc 4>a (tb))- The square of the field-field correlation function takes a simple form 
in the equivalent channel approximation: Assuming that (Tj a ) = (1/N 2 ) £\ (Tj a ) = g/N 2 , 
equation © factorizes and 

C 1 (A,1;B,2) = |F(r A ,r B )| 2 |F(r l7 r 2 )| 2 (7) 

where 2)| 2 can be written in terms of the Green function (j^J: 

In the large- iV limit, the Green function of the clean waveguide tends to the free-space Green 
function, exp(z/cr)/(47rr) in the case of a 3D conductor, and we have 

sin(/c|Ar ? i 2 |) 



\F(r 1 , r : 



|2 



fc|Ari 2 | 

and analogous results holds in 2D with the "sine" replaced by the Bessel function J 



(9) 



\F(r u r 2 )\ 2 ^ J 2 (k\Ar 12 \). (10) 

The behavior of the function |F| 2 versus kAr = k\xi — x 2 \ is illustrated in Figure 1. Black 
dots correspond to the large- N limit, eq. (jlOj) . For a finite width, the field correlation 
function (eq. EJ) strongly depends on the position of the detectors. Finite size effects can 
be minimized by using the experimental approach of ref. 3]: we consider the average of 
\F\ 2 when x\ and x 2 are uniformly distributed over the interval (W/10, W — W/10). As can 
be seen in Figure 1, the averaged \F\ 2 (dashed line) already presents the typical large- iV 
behavior for the case of N = 20. 
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FIG. 1: Square of the field correlation function for a single source as a function of the distance 
between detectors for a 2D waveguide. The black dots correspond to the expected large-iV limit 
(eq. I10JI assuming that all channels are equivalent. For finite N (N = 20), the results represent the 
average over different source and detector positions: The large circles are the results of microscopic 
numerical calculations (for g = 1.12) based on the model system sketched in the inset. The contin- 
uous line is the result of equations El and El with {T a j ) obtained from the numerical calculations 
while the dashed line corresponds to the equivalent channel approach, (T a j) = g/N 2 (eq. EJ). 

The spatial intensity correlation function can be denned as 

C(A, 1; B, 2) = (I (A, 1)I(B, 2)> - (I (A, 1)) (/(£?, 2)) (11) 
where the first term of the right hand side is given by 

(I(A, 1)I(B,2)> = E E {( c ^«') x 

aa'bb 1 ii'jj' 

(0;(ri)0+*(r 1 )0+(r 2 )0+*(r 2 )) (t ja t*, a ,t ib t*, b ,)} (12) 



In contrast with field-field correlations, the calculation of the averages presents subtle 
properties directly related to the symmetry and statistical properties of the S matrix. Calling 
Ti, T 2 , ■ ■ ■ i the set of N common transmission eigenvalues of the four Hermitian matrices 
tt^, tt', rr^ and ff^, the matrix S can also be written in terms of the T n 's by means of the 
polar decomposition lj, [l^ 





FIG. 2: Plots of the intensity correlation function. Inset: Long range intensity correlations (C2) as 
a function of the length s = L/i of the system. Black dots are the results of numerical calculations. 

where are unitary matrices and 

T = diag(7i, 72, ... , Tn) is a N x N diagonal matrix with the transmission eigenvalues on 
the diagonal. The transmission and reflection matrices can, then, be written as 



r ba = -XX" (v 7 ! - %) u 

n 

t- = V M (2) (Jt\vP> 



(13) 
(14) 
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One of the key assumptions in the macroscopic approach is the hypothesis of isotropy 



14, 



lq . Under this hypothesis the statistical distribution of the transmission eigenvalues 



{7~ n } is independent of the unitary matrices u^, and the calculation of the statistical averages 
in (|12|) factorizes. Moreover, u^ 1 ) and are statistically independent from each other, 
each being distributed according the invariant measure of the unitary group. By using 
the averages over the unitary group ((uj n )(uji n i)*) and ((uj n u im )(uj> n >Ui> m >)*) (evaluated by 
Mello in ref. [itI]). after some algebra, we find 



(15) 



(tjatfa'tibti'v) 



[A N (T 2 ) — B N (T 2 )] (Sij'Si'jSab'5 a 'b + Sjj'Su'Saa'Sbb') 

+ [A N (T 2 ) - B N (T 2 )] (Sjj'du^ab'da'b + 5ij>8 Vj 8 aa >8 hV ) (16) 
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introducing the notation g = (T) = ^2 n {%), 

var{g} = (T 2 ) — (T) 2 , (T 2 ) = ^ n {T 2 ) and the coefficients An and B^ given by 



A - N " + 1 1 - R 2 1 

^Jf — 77^ — 77T^ 77^7 i -^N 



AT2 (AT2_!)2 ' - jV(iV 2 -l) 2 

^From equation and taking i = = j', 
a = a',b = b', we easily recover the well known channel-channel correlation function Cj a i b 

{T ja T ib } 



(T 3a )(T tb ) 
where 



1 = CiO^) + C 2 (5 ab + 5 l3 ) + C 3 (17) 
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& = -— 2 [A N {T 2 ) - B N {T 2 )} (18) 

C 2 = —[A N (T 2 )-B N (T 2 )} (19) 
C 3 = d - 1 (20) 

The angular correlation function has the same structure as that found first by Feng et al. P| , 
the three coefficients C\, C 2 and C 3 corresponding, respectively, to short, long and infinite 
range correlations. 

In order to obtain the spatial correlation function, we substitute the expressions (|T5|) and 
(JTOJ) for the averages in equations (JOJ) and (|T2|) . 

(7(^1)7(6,2)) 
C(Ar a6 ,Ar u )= <J(flil))<J(6>2)) -l = 

d (|i^(r , r 6 )| 2 |i^(r l5 r 2 )| 2 ) 

+ C 2 {\F{7 ai r b )\ 2 + \F{r Xl T 2 )\ 2 ) + C z (21) 

where |F(1,2)| 2 is the square of the field correlation function (eq. IHJ) and the correlation 
coefficients Ci, C2 and C3 are exactly the same as those appearing in equation (|17|). The 
structure of the spatial correlations is, then, equivalent to that obtained for channel corre- 
lations with the angular u 5 ab " functions replaced by the spatial functions \F(a,b)\ 2 . 

The behavior of the spatial correlations is illustrated in Figure 2. For a single source and 
two detectors ( Ar ab = AR = and Ari 2 = Ar) , or viceversa, the correlation function (in a 
2D waveguide) is given by m (C\ + C 2 )(Jo(kAr)) 2 + C 2 + C3 and approach a constant as Ar 
increases. For large AT?, the correlations as a function of Ar behave as f=s C 2 ( Jn(k Ar)) 2 + 
C 3 . These results are in qualitative agreement with the results of reference (lOj. As a 
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FIG. 3: Correlation coefficients C\, C2 and C3 as a function of the length s = L/£ of the system. 
Continuous lines are the results of the DMPK equation. Dashed line represents the perturbative 
1/N expansion (after ref. [{]]). Circles are numerical results for the model system sketched in Fig. 
1 

matter of fact, our expressions f!21l) and (jHJ) ( or El in the large- N limit) are consistent, 
after a slight reordering of the terms, with the expression (5) in reference [lo| . which was 
obtained through a perturbative diagrammatic expansion. The equivalence with polarization 
correlations (equation (2) of reference Ql) is also evident. However, while the diagrammatic 
expansions are strictly valid in the diffusive regime, our result does not depend at all on 
the transport regime, and are a direct consequence of the isotropy hypothesis. Only the 
relative size of Ci, C2 and C3 will depend on the length of the system L and the mean free 
path £ through the distribution of transmission eigenvalues P({T n }, s), with s = L/£. Figure 
3 summarizes the dependence of the correlation coefficients with the leng th of the system 
as obtained from a Monte Carlo calculations of the DMPK equation 18j. In the diffusive 

nn 

regime (1 <^ g ^ N), the correlation factors take the well known values p, |6[: C\ ~ 1, 
gC2 ~ 2/3 and g 2 C% ~ 2/15. As the length of the system decreases the exact solution of the 
DMPK equation is well described by the known analytical results based on 1/N expansion 
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(dashed line in Fig. 3). While the correlations are positive in the diffusive regime, our 
results predict a transition to negative correlations for both the angular and the spatial case 
(see inset in Fig. 2). 

In order to confirm the predictions of the macroscopic approach, we have performed 
extensive numerical calculations based on the simple two-dimensional (2D) model sketched 
in Fig. 1. For electromagnetic waves, we assume s-polarization with the electric vector 



8 



parallel to the walls. The disordered region is divided in small rectangular regions of section 
5 Z x 5 X . Within each slice the refraction index n R has random values distributed uniformly 
in the interval (1 - 5n R , 1 + 5n R ). We take 5 X = W/10, 5 Z = W/20 and 5n R = 0.025, for a 
20-mode waveguide (W/X = 10.25) as in Figure 1. Transmission and reflection coefficients 
are exactly calculated by solving the 2D wave equation by mode matching at each 5^-slice, 
together with a generalized scattering-matrix technique |l9j]. Figure 1 shows the behavior 
of the averaged field correlations (for s = 15, (g) = 1.12). The numerical results show that 
transport is not fully isotropic (the distribution of the transmission coefficients Tj a has some 
dependence on the particular values of j and a). The numerical intensity correlations are 
very close to the expected behavior (see Fig. 2) although they are not fully described by 
equations H7I and I2T1 However we can extract the correlation coefficients C\,Ci and C3 
from a least square fitting of spatial or angular numerical correlations with equation (|2~Tj) 
or (|17|) respectively. In both cases, the obtained coefficients (circles in Figure 3) are in full 
agreement with the predictions of the DMPK approach. 
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the European Integrated Project "Molecular Imaging" (LSHG-CT-2003-503259). One of us 
(GC) wants to acknowledge the hospitality of the D.F.M.C. (U.A.M. - Spain) where this 
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